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Abstract 

We propose a simple dynamical model that generates networks with power-law degree distri- 

c3: 

. butions with the exponent 2 through rewiring only. At each time step, two nodes, i and j, are 

randomly selected, and one incoming link to i is redirected to j with the rewiring probability R, 
determined only by degrees of two nodes, ki and kj, while giving preference to high-degree nodes. 
To take the structure of networks into account, we also consider what types of networks are of 

o , 

O , interest, whether links are directed or not, and how we choose a rewiring link out of all incoming 
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links to i, as a result, specifying 24 different cases of the model. We then observe numerically that 

> ' 

networks will evolve to steady states with power-law degree distributions when parameters of the 
I/"") , model satisfy certain conditions. 
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Complex systems have been studied in various research fields as a general framework to 
describe systems that usually consist of interacting elements, and, in most cases, these sys- 
tems as a whole exhibit emerging behavior. We can represent complex systems as networks 
with nodes (or vertices) and links (or edges), where nodes represent elements of a given 
system, and links represent interactions (physical or conceptual) between two elements. 
Recently, structural and dynamical properties of these networks and their underlying mech- 
anisms have become the focal point of many researches Q^E,!^. One well-known property 
of many natural and artificial networks alike is the power-law behavior of degree distribu- 
tions P(k) (the degree k of a node is the number of connected links), and such networks are 
called scale-free networks; in other words, P{k) oc A; -7 where 7 is the exponent, usually in 
;he range of 2 < 7 < 3. Examples are citations of scientific papers |5|, links in web pages 
f|, the Internet 3|, metabolic networks ^, and brain functional networks jj, and so on. 

Many generating mechanisms have been proposed so far to explain this phenomenon. One 
popular approach is using stochastic dynamical processes, in which networks evolve with time 
and approach scale-free states, ensembles of network structures sharing specific statistical 
properties. If we classify these mechanisms into two classes, one is growing mechanisms, 
where both the number of nodes, N, and the number of links, E, increase with time, and 
the other is non-growing mechanisms, where N and E are either constant or varying but 
limited to finite values. One of the most studied growing mechanisms is the preferential- 
attachment (PA) mechanism , where a node and a fixed number of links connected to this 
node are added at each time step, while other ends of these links are chosen from existing 
nodes with the probability proportional to their degrees. There are many generalizations 
and extensions of PA mechanism, usually with additional features: for example, rewiring 



12j, linking between existing nodes [13J], age and deactivation 



11 1 , attractiveness of nodes 

14 1, local knowledge Q], and link weig hts [16J. While PA mechanism is valid to many 
scale-free growing networks, it is not suitable for non-growing networks, and there also are 
numerous non-growing mechanisms that use some additional features or information: for 
e, transitive linking and aging 17], fitness of nodes Q|, adding links using nonlinear 
19], local optimization using memories of nodes j^j], optimization of Hamiltonians [21], 
and merging and regeneration of nodes (2^] . Models mentioned so far share the property of 
power-law degree distributions, but they also attempt to capture different aspects of various 
types of real networks, which we don't describe details here. 
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In this paper, we propose a non-growing mechanism (both N and E are constant) with a 
stochastic preferential-rewirinig process, where only degree information of interacting nodes 
is used for the dynamics. In other words, nodes and links have no variables that represent 
memory, attractiveness, or weight, and the global structural knowledge of the network is not 
needed, thereby making the dynamics simple and local. Our network model is an application 
of the general model introduced in Ref. |23j|, which we first describe briefly. 

The general model deals with iV-element systems, where each element i is represented 
by a quantity ki, a non- negative integer. At each time step, two elements, % and j, are 
randomly chosen out of iV elements and one unit of quantity is given from i to j with the 
exchange probability R, determined only by ki and kj. When the exchange occurs, ki and 
kj are changed to ki — 1 and kj + 1 (the average quantity, a = J2iLi h/N = (k), is constant). 
The exchange probability is defined using a parameter (3 as below, 



R 



1 (0 < ki < kj) 

(3 ih>k j ) , (i) 

(ki = 0) 



where <J3 < 1. This is a generalized zero-range process (ZRP) in a fully connected 
geometry 2^, even though R as a function of fc, and kj in Eq. (JIJ is rather unique. It was 
found analytically and numerically that the distribution of fcj's exhibits the power law when 
a and j3 satisfy the condition in the limit of large N, 

«c(/5) = r 4^1n[e- 1 ^] <a, (2) 

where e is a small positive constant, estimated to be about 10 -3 . This equation defines a 
power-law regime inside (a, /3)-space; i.e., to have a power-law distribution when (3 is given, 
a should be greater than the critical value, a c ((3) [or when a is given, (3 should be less than 
the critical value, (3 c (a), which is the inverse of a = a c ((3)]. For example, from a c (0.5) ~ 6, 
we deduce that (3 should be less than 0.5 to have a power-law distribution when a = 6. 
When Eq. (J2J) is satisfied, the stationary probability distribution, P(k), and the stationary 
cumulative distribution, P(>k), have been found analytically as below, 

(3 1 



P(k) 



1 -fi [k + 1/(1-/3)} [k + (3/(1- (3)Y 



FIG. 1: A diagram describing the rewiring dynamics. Links of the given network can be either 
(a) undirected, or (b) directed. For both cases, two nodes, i and j, are randomly chosen (ki = 2 
and kj = 3), and the node i' is one of the neighbors of i. (For networks with directed links, fcj 
represents the in-degree of the network, and %' is one of the neighbors with an outgoing link to i.) 
Rewiring from (i 1 , i) to (i',j) can occur with the rewiring probability R. 

which is the power-law distribution with 7 = 2 when k(l — (3) ^> 1. 

A network with undirected links (a = 2E/N) is an example of such A^-element systems, 
where ki represents the degree of a node % and R is the rewiring probability (for networks 
with directed links, ki represents the in-degree, the number of incoming links, of a node and 
a = E/N). In Fig. ^ the rewiring dynamics is described using a schematic diagram. A 
rewiring process involves three nodes: % and j are randomly chosen, and i', a pivot node, 
is a neighbor of i, which means a node directly connected to i. When ki and kj are given, 
rewiring from to can occur with the rewiring probability R in Eq. (JTJ. Due to 

the presence of the pivot node and different network types, our network model should be 
more than just a ZRP if we consider the link structure of networks, and we introduce three 
additional distinctions: how to choose %' among neighbors of % (pivot type), whether links 
are directed or not (link type), and whether loops and/or multiple links are allowed or not 
(network type). (When two end nodes of a link are the same, we call it a loop, and when 
there are more than one link between two nodes, we call them multiple links). 

Here we introduce three pivot types: (i) the node with the smallest degree ('5"), (ii) the 
randomly chosen node ('i?'), (iii) the node with the largest degree ('£'). There are two link 
types: (i) undirected ('£/'), (ii) directed ('-D'). And there can be four types of networks: (i) 
both loops and multiple links are allowed ('1'), (ii) only loops are allowed ('2'), (iii) only 
multiple links are allowed ('3'), (iv) neither a loop nor multiple links are allowed ('4'). Then 
there are 24 different cases of the model. For convenience, we will use a three-letter notation 
from here on, where the first letter represents the pivot type, the second letter represents 
the link type, and the last letter represents the network type. For example, SD4 represents 
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FIG. 2: Stationary cumulative degree distributions, P(> k), of 12 different cases, when a = 6, 
/? = 0.5, and iV = 10000: (a) SU1 and SUA, (b) 5D1 and ST>4, (c) RU1 and AZ74, (d) RD1 and 
RD4, (e) LC/1 and L{74, (f) L.D1 and LD4. The thin solid line represents the analytic solution, 
P(>k) = 1/(1 + k), of the general model when /? = 0.5 as iV — ► oo. 



networks with directed links, which allow neither a loop nor multiple links, when i' is the 
neighbor of % with the smallest degree. In addition, three more details should be mentioned: 
(i) while choosing the pivot node, if there are more than one node with the smallest or largest 
degree, one node is chosen randomly out of them, (ii) when the link type is D (directed), 
out-degrees of neighbors are used when choosing the pivot node, (iii) when the network type 
is not 1, a rewiring attempt can be aborted if it results in a loop or multiple links. 

In Fig. El we observe stationary degree distributions, P(>k), of 12 different cases when 
a = 6, (3 = 0.5, and iV = 10000. If we look at the cases when both loops and multiple links 
are allowed (network type 1), degree distributions are the same as those from the general 
model, because degree distributions don't depend on the pivot type and the link type even 
though other structural properties of the network do. These cases can be regarded as ZRPs 
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when we are only interested in degree distributions. On the other hand, when neither loops 
nor multiple links are allowed (network type 4), the degree distributions are different from 
the general model because rewiring attempts can be aborted depending on the current link 
structure of the network; as a result, cutoffs are usually lower. When the pivot type is S, 
the distribution is closer to that of the general model even though cutoffs are still somewhat 
lower, because the probability of aborting a rewiring attempt is lower when the pivot node 
has the smallest degree [see Figs. Efa) and (b)], and that is especially true for cases with 
undirected links. When the network type is 2 or 3 (not shown here), results are almost the 
same as those of the network type 4 or 1, respectively. This shows us that the presence of 
loops does not have much impact on stationary distributions. 

In Fig. El we observe the behavior of the distributions of the SUA and SDA networks 
more closely. We vary one parameter while the other two parameters are fixed. In Figs. 
Ota) and (b), a is varied when N and (3 are fixed, while in Figs. |3fc) and (d), (3 is varied 
when N and a are fixed. We observe that the stationary degree distributions are close to 
the theoretical curve when a is close to a c {(3) [or (3 is close to f3 c (a)] unlike the general 
model, where stationary distributions were found to follow the power law when a > a c (f3) 
[or (3 < (3 c {a)]. In Figs. Ofe) and (f), N is varied when a and (3 are fixed. We observe that 
cutoffs increase as N increases, because the sparser the network gets, the smaller the chance 
of aborting the rewiring attempt becomes. 

In Fig. |H we observe other structural properties of the SU 4 network of the steady state 
when a = 6, (3 = 0.5, and iV = 20000. There are several well-known statistical measures 



such as the clustering coefficient, C 12511 . the average shortest path length, L, and the 
Pearson correlation coefficient of degrees, r |26j . Here an initial network is the Erdos-Renyi 
(ER) random network, where C = 0.0003, L = 5.73, and r = 0.00008. When the stationary 
distribution is reached at T/N ~ 10 5 , we find that C = 0.04(0.002), L = 4.53(0.03), and 
r = 0.2(0.02), showing the small-world behavior [27]. These networks have more than one 
component: for ER networks, the size of the largest component is 19958 and P(0) = 0.002, 
and they become 10199(53) and 0.467(0.003) respectively in the steady state. In Fig. Ufa), 
the average clustering coefficient of nodes with degree k, C(k), is found, and we can observe 
that C(k) increases with k. And the inset shows how C increases with time until T/N ~ 10 5 . 
In Fig.H^b), the averaged degree of nearest neighbors of nodes with degree k, (k nn ), is found, 
and we can observe that (k nn ) increases with k also. The inset shows how r changes with 
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FIG. 3: The stationary cumulative degree distributions, P{>k), of networks, SUA and SD4, when 
one of three parameters of the model, a, (3 and N, is varied, (a), (b) a is varied; (c), (d) (3 is varied; 
and (e), (f) N is varied. The thin solid line represents the analytic solution, P{>k) = 1/(1 + k), 
of the general model when (3 = 0.5 as N 
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time, and both results imply that this network is assortatively mixed. In both cases, the 
increasing behavior stops at k ~ 200, and this is probably due to the finite-size effect [the 
distribution deviates from the theoretical curve at k ~ 200 in Fig. Ete)]. 

In conclusion, we introduced the rewiring model that generates scale-free networks with 
7 = 2. We showed numerical results for 12 different cases of the model, and observed that, 
when the network can have loops and multiple links, stationary degree distributions for all 
cases, including those from the general model, follow the same power law if the condition 
a > a c {(3) is satisfied. While there exist networks with multiple links (for example, citation 
network), networks without a loop and multiple links have been of interest in most studies. 
We also observed that, when the network can have neither a loop nor multiple links, scale- 
free networks are generated if the condition a ~ a c (j3) is satisfied. This model was not 
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FIG. 4: Structural properties of stationary SUA networks when a = 6, (3 = 0.5, and N = 20000. 
Data points (•) are obtained from averaging values of individual nodes in log-sized bins in both 
cases, (a) The averaged clustering coefficient, C(k), versus A; in a log-log plot at T/N = 10 6 . Inset: 
The clustering coefficient, C, versus T/N in a log-log plot, (b) The averaged degree of nearest 
neighbors of nodes with deg P66 (fenn) •) V6rSUS k in 8l log-log plot at T/N = 10 6 . It shows that 
the stationary network is assortatively mixed. Inset: The Pearson correlation coefficient of the 
degrees, r, versus T/N in a semi-log plot. 



conceived to describe a certain type of real networks, but it was from an effort to find a 
simple rewiring model generating scale-free networks, with intention of using local degree 
information only and incorporating the rich-get-richer phenomenon. As one of the simplest 
rewiring models for scale-free networks, this model can be extended to meet specific needs 
of some real non-growing networks. 
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